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ABSTRACT 

When present, extended disks of neutral hydrogen around spiral galaxies show a 
remarkably uniform velocity dispersion of ~ 6 km s~^. Since stellar winds and supernovae 
are largely absent in such regions, neither the magnitude nor the constancy of this 
number can be accounted for in the classical picture in which interstellar turbulence is 
driven by stellar energy sources. Here we suggest that magnetic fields with strengths 
of a few microgauss in these extended disks allow energy to be extracted from galactic 
differential rotation through MHD driven turbulence. The magnitude and constancy of 
the observed velocity dispersion may be understood if its value is Alfvenic. Moreover, 
by providing a simple explanation for a lower bound to the gaseous velocity fluctuations, 
MHD processes may account for the sharp outer edge to star formation in galaxy disks. 

Subject headings: galaxies: ISM — galaxies: kinematics and dynamics — hy- 
drodynamics — instabilities — radio lines: galaxies — turbulence 
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1. Introduction 

It has long been recognized that the intersteUar 
medium (ISM) in galactic disks is turbulent (e.g. 
Scalo 1987, Dickey & Lockman 1990). Within our 
own Galaxy, line widths of individual molecular and 
H I clouds, and cloud complexes greatly exceed the 
expected thermal width, when kinetic temperatures 
can be estimated (Miesch & Bally 1994; Caselli & 
Myers 1995). There is also evidence for a high ve- 
locity dispersion population of H I clouds in both 
the Milky Way Galaxy (Radhakrishnan & Srinivasan 
1980, Kulkarni & Fich 1985), and for a high veloc- 
ity tail in the gas where the disk is optically bright 
in some external galaxies (Dickey, Hanson, & Helou 
1990; Boulanger & Viallefond 1992; Kamphuis & San- 
cisi 1993; Schulman et al. 1996; Braun 1997). 

The classical picture (e.g. Spitzer 1978) is that su- 
pernovae and other stellar processes (e.g. winds, out- 
flows, etc.) supply the requisite energy to maintain 
turbulent cloud velocities against dissipative losses. 
In the absence of such sources, one would therefore 
expect interstellar turbulence to decay rapidly. In 
particular, turbulent motion should be negligible be- 
yond the outer optical edges of spiral galax;ies, where 
there are few stars. 

This expectation is not in accord with observa- 
tions. Radio observers report that 21 cm dispersions 
do not drop below 5-7 km s~^, regardless of how low 
the optical surface brightness becomes (e.g. Dickey et 
al. 1990; Kamphuis 1993). If this line width is a mea- 
sure of the turbulent motion in an ensemble of cool 
clouds, the medium is being stirred by other means. 
We here propose that MHD instabilities (Balbus & 
Hawley 1991) in the differentially rotating gas layer 
are responsible for a minimum level of turbulence in 
all gas disks. 

The need for MHD mediation of interstellar tur- 
bulence has been noted by a number of other inves- 
tigators. Magnetic fields in galactic disks tend to be 
partly tangled on small scales, where lines of force 
have been observed to thread and be compressed by 
dense clouds, but also to retain some large-scale co- 
herence (Heiles et al. 1993; Heiles 1996). On the 
largest Galactic scales, the field energy density is less 
than that of the random cloud motions, but proba- 
bly in excess of the average thermal energy density. 
Within molecular cloud complexes, the field strength 
can be dominant, and MHD processes in these sys- 
tems are central to understanding their dynamics and 



evolution (e.g. McKee & Zweibel 1995; Gammie & 
Ostriker 1996). 

The extended H I disks of some galaxies are a rel- 
atively clean laboratory in which interstellar turbu- 
lence can be observed on its largest scales. It is clearly 
present in these systems, even though star formation 
and the associated supernova rates are very low. If 
the turbulent motions in this regime are Alfvenic, 
we show that they are maintained through magnet- 
ically mediated dynamical heating from differential 
rotation. We present a derivation of this heating 
rate when magnetic fields and self-gravity are both 
present, though in our specific application the pro- 
cess is magnetically dominated. Our result is applied 
to the particular case of NGC 1058, a face-on disk 
galaxy with a well-studied extended H I disk. Both 
the constancy and magnitude of the observed velocity 
dispersion follow naturally from an MHD explanation. 

2. Turbulence in extended H I disks 

A compilation of high-quality neutral hydrogen line 
width data on several galaxies may be found in Kam- 
phuis (1993). The principal conclusions from this 
study are that 21 cm dispersions arc cr ^ 10 — 12 km 
s~^ in the bright optical disk (with evidence for still 
higher values in spiral arms), and that the dispersions 
decline to no less than 5 — 7 km when the opti- 
cal surface brightness falls below 25 B mag arcsec"^. 
The higher spatial resolution data presented by Braun 
(1997) indicate that broad lines are foimd in a "high 
brightness network" absent in the optically faint outer 
disks. 

Dickey et al. (1990) drew particular attention to 
the remarkable lack of spatial variation of the H I 
line width in the gas layer beyond the optical disk 
in NGC 1058. Data from other galaxies show a sim- 
ilar trend, and indicate that the velocity dispersion 
puzzle is widespread. Both NGC 1058 and NGC 
5474 (Rownd et al. 1994) are almost face-on, allowing 
clean measurements of line widths with little contri- 
bution from beam smearing. The velocity dispersions 
of these galaxies appear to approach a steady value of 
cr ~ 6 km s~^ at radii outside the bright disk. (The 
missing flux in their observations could, in principle, 
affect this result, but a similar dispersion is obtained 
for other galaxies (e.g. Kamphuis 1993) in which the 
full flux is recovered.) Dickey et al. stress that where 
the line width is small, the H I line profiles are almost 
perfect Gaussians. 
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Let us focus now upon NGC 1058, a small galaxy 
with an adopted distance of 10 Mpc. Its optical ra- 
dius R25 = 90" or 4.5 kpc and the H I layer extends 
to 220" or 11 kpc. Since it is almost facc-on, its ro- 
tation curve is poorly known; we adopt a flat circular 
velocity of 150 km s~^, consistent with the H I veloc- 
ity field with 0.05 ^ sini ^0.1. Dickey et al. report 
a H I column density ~ 3 x 10^" cm~^ at i? ~ 5 kpc 
falUng to 2 X 10^° cm'^ at i? ~ 10 kpc. These are 
likely to be underestimates, however, since their intcr- 
ferometric data miss about half the flux determined 
from single dish observations. Therefore, as fiducial 
values, wc double the column density numbers in our 
calculations below, adjust for helium content and also 
assume there is no significant quantity of molecular 
gas in the outer layer. We adopt a total thickness 
of the outer H I layer of h = 400 pc - about twice 
the thickness of the gaseous disk of the Milky Way 
in the Solar neighborhood, since wc expect the layer 
to flare to some extent. We summarise the properties 
we adopt in Table 1. 

The constant internal H I line width, Dickey et al. 
argue, could represent either the thermal temperature 
of the gas or the turbulent motion of many small, 
cooler clouds. Neither interpretation is free of puzzles. 

If the observed 1-D dispersion of ~ 5 — 7 km 
s~^ is thermal, it would imply a gas temperature of 

3 000 — 6 000 K, perhaps less if there are some con- 
tributions to this width from beam smearing (small 
in a nearly face-on galaxy) and turbulence. Dickey 
(1996) points out that atomic hydrogen at a temper- 
ature of a few thousand degrees would be in an un- 
stable state. Braun (1998) expects the low ambient 
pressure to cause all gas to be in the warm phase at a 
temperature of ~ 10 000 K (cr ~ 9 km s""'^), which is 
inconsistent with the observed dispersion. While the 
principal heating and cooling processes at work in the 
Solar neighborhood (Wolfire et al. 1995) arc likely to 
be different in the far outer disk, they would have to 
be drastically so for gas to be in a stable phase at 
^ 5 000 K. Moreover, stars are forming at a very low 
rate in at least some of these gas layers (Ferguson et 
al. 1998), which must imply the existence of cooler, 
denser gas from which they can form. 

We therefore prefer to attribute the observed line 
width largely to turbulent motion in an ensemble 
of small, cool clouds, which immediately raises two 
related questions: i) What process accelerates the 
clouds, maintaining velocities well in excess of the in- 
ternal sound speed of individual clouds?; and ii) Why 



should the turbulent velocities be so uniform? We 
offer answers to both these questions in §4. 

2.1. Supernova heating 

Supernovac (SNc) are generally assumed to be the 
dominant energy source for the interstellar turbulent 
cascade (e.g. Spitzer 1978; Norman & Ferrara 1996). 
The stellar density and supernova rate (SNR) in the 
gas layer beyond the Holmberg radius is clearly much 
lower than in the bright inner parts of galaxies, but 
could it still be sufficient to drive the turbulence? 

Deep CCD images of NGC 1058 (Ferguson et al. 
1998) have confirmed the very faint outer spiral arms, 
first detected by Tamman (1974) in his search for the 
stellar population responsible for the two SNe de- 
tected during the 1960s in the outer parts of this 
galaxy. Ferguson (1997) estimates an azimuthally 
averaged surface brightness declining to 28 B mag 
arcsec"^, with a regular spiral pattern. Ferguson et 
al. also detect small H 11 regions along these spiral 
arms, but the inter-arm region seems to be almost 
devoid of any star formation. 

From the Ha flux, Ferguson (private communica- 
tion) estimates the azimuthally averaged star forma- 
tion rate (SFR) per unit area at ^ 1.5 times the op- 
tical radius to be ~ 5 x 10~^^ Mq pc^^ yr^^. A high 
estimate for the SNR for this SFR is lO^i^ SN pc'^ 
yr~^ (e.g. Leitherer & Heckman 1995). Assuming me- 
chanical energy input of elO'^^ ergs per SN, where e 
is the highly uncertain efficiency factor, we obtain an 
energy input rate of~ e3 x 10^^''(400pc//i) ergs cm^'' 

For e ^ 0.01 (Chevalier 1998, private commimica- 
tion), this energy input rate, while highly uncertain, 
is not decisively less than that of the differential rota- 
tion source we consider below. Neverthless, it seems 
unlikely to be the principal source of cloud motions, 
because there is no correlation whatsoever between 
the SFR and H I velocity dispersions across the ex- 
tended disk. Dickey et al. (1990) stress the uniformity 
of the dispersion, while Ferguson et al. observe H 11 
regions to lie almost exclusively in very narrow arms. 

Braun's (1997) high spatial resolution studies of 
H I in the bright inner parts of other galaxies show a 
"high brightness network" of distinctly non-Gaussian 
line profiles with broad, high velocity tails. He further 
reports (Braun 1998) a loose correspondence between 
the network of extended H 1 linewidths in M31 and 
the observed Ha emission, and argues for a causal 
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connection between linewicltli and energy injection 
from stellar activity. Thus turbulence driven by stel- 
lar activity appears to create broad wings in the line 
profiles in localized regions where the energy is de- 
posited. It is hard to see how this activity could also 
produce the observed uniform level of near Gaussian 
line profiles on larger scales. 

2.2. Other sources of turbulence 

Other sources of turbulence can be imagined. In- 

fall models include the returns from a galactic foun- 
tain (Bregman 1980; Schulman et al. 1996), or chim- 
ney (Norman & Ikeuchi 1989), or simply the direct 
accretion of external intergalaetic matter (Toth & 
Ostriker 1992; Kamphuis & Sancisi 1993). Gravita- 
tional scattering by transient spiral waves (Carlberg 
& Sellwood 1985; Jenkins & Binney 1990; Toomre & 
Kalnajs 1991) is yet another possible energy source. 

Dickey et al. (1990) already remarked that the very 
low SFR in NGC 1058 is almost certainly inadequate 
to drive a vigorous galactic fountain. Even granting 
the presence of a weak fountain, the measured H I 
dispersion is no lower than in the outer parts of other 
galaxies with considerably more active star formation 
(Kamphuis 1993). Furthermore, stirring of the H I 
layer by infalling dwarf galaxies and debris is unlikely 
to maintain a uniform level of mild turbulence every- 
where, and should result in radiative emission of a 
substantial fraction of the infall energy. 

Transient spiral waves produce a choppy potential 
sea which scatters stars, and any other material in the 
disk, causing random motion to rise. Where most of 
the disk mass is in the collisionless component - i.e. 
the stars the process is self-limiting, since the spiral 
waves quickly become weaker as the velocity disper- 
sion of the stars rises. Because gas is able to dissipate 
turbulent energy through inelastic collisions, its fate 
will be different, however. It seems unlikely that gas 
would settle to a smooth distribution on large scales 
with a level of turbulence resulting from a balance 
between scattering and dissipation; it is more likely 
that stellar spiral arms promote the formation of large 
gas concentrations in which further gravitational in- 
stability will lead to star formation. In parts of the 
disk in which gas is the dominant mass component, if 
gravitational instability is present on large scales, we 
might expect it to cascade directly to forming stars. 

The importance of self-gravity in a disk may be 
determined by evaluating the usual local stability pa- 



rameter (e.g. Binney & Tremaine 1987) 



where thickness corrections have been neglected. Here 
a is the velocity dispersion, the epicyclic frequency 
K = \/2Vcitc/ R for a flat rotation curve, and E is the 
surface density. Values of this parameter for the ax- 
ially symmetrized gas only in the outer H I layer of 
NGC 1058 arc evaluated in Table 1. The Q value 
for the extremely faint stellar disk seems to be huge: 
Q ~ 60(6kpc/i?)(a„/10kms"^)/TB, where is 
the mass to B-band luminosity in Solar units; any 
supporting response from the stars must be utterly 
negligible in this region. 

The values of Q given in Table 1 are uncertain and 
rather modest to argue strongly that the gas layer is 
clearly stable. We note, however, that the gas sur- 
face density varies by perhaps a factor two between 
the peaks in the spiral arms and the inter-arm level. 
Ferguson et al. (1998) report mild star formation, a 
clear indicator of local Jeans instability, in the arms 
but the paucity of detectable star formation between 
the arms suggests that the layer is stable there. The 
origin of the spiral arms in this outer layer is unclear, 
but gravitational instability seems quite unlikely since 
the Jeans length (Acrit = 47r^GSgas/K^) is much too 
small for such a large-scale and symmetric two-armed 
spiral (Table 1). We conclude that gravitational in- 
stabilities are unlikely to be stirring the gas layer to 
maintain the low level of turbulence in the gas. 

3. Turbulent heating by coupling to differen- 
tial rotation 

In this section we consider the local turbulent dy- 
namics of a galactic disk. Disk material orbits in the 
global potential of the galaxy, but is locally subject 
both to a magnetic field and (in principle) to its own 
self-gravity. We require the form of the turbulent vol- 
umetric heating rate under these conditions. For the 
reasons just given, self-gravity is not of direct impor- 
tance for the gaseous disks under consideration here 
but we include it in the analysis of this section for 
the sake of completeness, and because it is an inter- 
esting problem in its own right. In the presence of 
self-gravity, the volumetric turbulent heating rate of 
the gas due to differential rotation is a direct gener- 
alization of the nonself-gravitating case: a coupling 
between the effective stress tensor and the large scale 
angular velocity gradient. 
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In a homogeneous gas, if the field is subthemial, 
the free energy of differential rotation drives a dy- 
namical instability (Balbus & Hawley 1991, 1998), 
extracting rotational energy and depositing it in tur- 
bulent motions. The physical mechanism for the gen- 
eration of turbulent motion is clear. Mass elements 
orbiting in the fluid at slightly differing radii, but cou- 
pled magnetically, pull on each other as the shear at- 
tempts to separate them as if connected by a weak 
spring. The effect of their mutual forces is to remove 
angular momentum from the one which has less and 
donate it to the one which has more. As a result, 
the radial separation of the elements increases, which 
increases the difference in angular velocity and the in- 
stability runs away - provided the spring is not strong 
enough to resist. We have a rather different case in 
mind here: a highly inhomogeneous gas in which the 
relevant "thermal" motions refer to the macroscopic 
velocity dispersion. It is likely that random magnetic 
stresses will continue to tap into the differential ro- 
tation as a source of turbulence for relatively strong 
fields as well (e.g. Eardley & Lightman 1975). Indeed, 
the expected outcome of differential rotation and any 
radial magnetic field component is a positive radial- 
azimuthal Maxwell stress; this alone is sufficient to 
drive noncircular velocity ffuctuations (cf. below). 

Adopting a standard (i?, 4>, z) cylindrical coordi- 
nate system, we denote the circular velocity as RQe^, 
where is a unit vector in the azimuthal direction. 
The velocity vector u is the difference between the 
true velocity and the azimuthal circular velocity; it is 
a fluctuation velocity satisfying u <C i?r2. We allow 
for the possibility of a mean, slowly varying drift ve- 
locity, denoted (u) , which is much less than the RMS 
fluctuation (it^)^/^. The angle brackets ( ) denote 
local averages in radius and height, but a complete 
average in azimuth. Thus 



|(u)| « {uY^^<t:Rn. 



(2) 



The Alfven velocity associated with the magnetic field 
B is given by 

"A = -?= (3) 



where p is the mass density. We do not make a 
formal distinction between the mean magnetic field 
and its fiuctuating component, but like its kinetic 
counterpart, the magnitude ua is assumed to satisfy 
UA < R^- 

We work in the standard local approximation in 
which R is assumed large enough that we can ignore 



geometrical curvature terms. Our starting point is the 
energy equation for the u fiuctuations with the grav- 
itational contribution written explicitly as a power 
term (see equation [89] in the review of Balbus & 
Hawley 1998): 



u\)) 



V-() = 



d\nR 
- {pu ■ V$) + 

'^/pv\Vu, 



UABUA4,)) 

PV • u) 



An 



(4) 



Here P is the gas pressure, v the kinematic viscosity, 
77 the resistivity and the operand in the divergence 
term V • ( ) is the energy fiux 



B 



(ipw^ + P)uH x(uxB 



(5) 



The gravitational potential <I> includes both a contri- 
bution from self-gravity, plus a contribution from the 
external large scale potential. Locally, the latter may 
be taken to depend upon z only, 

<i> = $sg(P, </>, z) (self gravity) -I- ^{z)cx (external). 

(6) 

Equation (Q) balances changes in the energy den- 
sity against the net fiux, and explicit sources and 
sinks. The first term on the right side is the energy 
released by the Reynolds-Maxwell stress, and is the 
key term in this paper. The other terms are the power 
extracted from the gravitational potential, the work 
done by pressure forces, and the viscous and resistive 
losses. Internal sources for interstellar energy fluctu- 
ations (e.g. SNe) can be included in the PV • u source 
term. 

We now show that the gravitational term can be 
manipulated into the form of a coupling to the stress 
tensor. Start with 



pu • V$ = V • (pu$) - $V • (pu). 



(7) 



Mass conservation and the Poisson equation then im- 
ply, 

-.fV.,pu,.*(_ + 0_)_ (8) 

(N.B. The second term in the brackets arises because 
the total velocity is u -f Pfie^.) Interchanging the 
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order of partial derivatives and integrating by parts 
leads to 



d V2$ 1 



dt A-kG AttG 



_ , ^av$\ 1 d ,„^,2 

V • $ V$ 

' dt J 2 9i ' ' 



(9) 



Similar, but more lengthy, manipulations lead to 



d 



1 



47rG 47rG 



1^^/ Ulni? 

--(•••) 



(10) 



The final partial (p derivative will vanish upon aver- 
aging, and is not explicitly written. Carrying this 
average through and returning to equation leads 
to 



{pu ■ V<I>) 



d /|V$P 



dt \ SttG 
'd 



-V • ( pu$ + $ ( — + 



47rG 



+ {PUGRUG4>)- 



d\nR 

where the gravitational velocity ug is defined by 



(11) 



Ug = 



(12) 



If we now place our findings back into the energy 
equation (^), we may write the result as 



dt 



dfl 
d\nR 



R4 



where the stress tensor Tr^ is defined as 

Tr^Ii = {p{uRU^ - UARUAcji + UGRUG4,)) 

and the supressed energy flux has become 



1 \ B 

-pu^ + P|uH x(uxB)+ pu$ 

2 J An 



(14) 



(15) 



The utility of writing the self-gravity in the form 
of a stress tensor was first noted by Lynden-Bell & 



Kalnajs (1972). We see that the most general stress 
tensor is a simple sum of the Reynolds, Maxwell, and 
Newtonian stresses. Free energy can be extracted 
from differential rotation at the rate per unit volume 
TR^dn/d\n R. 

It should be noted that the instabilities, both MHD 
and gravitational, drive the in-plane components of 
the velocity dispersion and do not couple directly to 
motions normal to the plane. Some degree of velocity 
anisotropy is therefore expected in turbulence driven 
by these stresses. 

4. Application to NGC 1058 

Differential rotation supplies energy to fluctuations 
at a rate given by 



R4 



dVt 
d\nR 



R(t> 



(16) 



for a flat rotation curve. We do not know Tr^ a pri- 
ori, but numerical simulations of magnetized nonself- 
gravitating disks, carried out under a wide variety 
of field geometries, equations of state, and numerical 
grids consistently yield 



R4 



0.6 



21 

8n 



(17) 



(Hawley, Gammie, & Balbus 1995), where the mag- 
netic field is evaluated at the time of saturation. 
A somewhat delicate point is how to relate to 
the mean magnetic field. The simulations reveal a 
field dominated by its largest scales, and we shall as- 
sume that mean field and fluctuations are compara- 
ble. Assuming the outer disk in NGC 1058 is not self- 
gravitating (as argued in §2.2), equation ( p7|) gives a 
simple, physically sensible estimate for the stress ten- 
sor, which we shall adopt. 

The scale for the fluctuating kinetic velocities is set 
by the Alfven velocity: 

{puRU^) ~ 0.6— + {puARUA^) (18) 

OTT 

Setting the left hand side to a fraction /i of pa^, 
where p is the mean density and the measured one- 
dimensional velocity disperion, and the entire right 
hand side to a fraction /2 of B /8tt, where B is the 
mean magnetic field, we find for a velocity dispersion 
of 6 km s~^ 



B = 3 /_tG 



/i 



1/2 



P 



10-24 g cm-3 



1/2 



(19) 
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The implied field strength is a few microgauss - a 
reasonable number for a spiral galaxy (e.g. Beck et 
al. 1996). Observational confirmation of this value 
would be challenging. [The radio continuum emission 
from NGC 1058 is very weak (van der Kruit & Shostak 
1984) because of its very low SFR.] 

If we turn the problem around and ask how a 
should vary with B and p, as the density varies be- 
tween and within the spiral arms say, a natural expla- 
nation for its near constancy emerges. Large scale but 
otherwise random field lines tend to follow a i? oc p" 
scaling, by making the usual assumption of flux freez- 
ing. The index n = 2/3 for isotropic spherical com- 
pression, while n = \ for compression in a plane. 
Since a is of order the Alfven speed, this implies 

a~(p)"-^/^ (20) 

i.e. a rather weak dependence on density. Therefore, 
a magnetic basis for the velocity dispersions gives, in 
addition to a sensible inferred field strength, a very 
simple basis for understanding the near constancy of 
(7 in gas-dominated galactic disks. 

Finally, let us estimate the energy deposition rate, 
using the numerical result (|l7|), and compare it with 
the SNR from §2.1. Inserting the above derived field 
strength of 3 pG, we obtain 

Tr^Q. ~ 1.4 X 10-2* ergs cm-^ s"^ (21) 

at i? = 7 kpc. Thus differential heating is certainly 
competitive with the 1% efficient SNR of §2.1, and 
provides a more natural explanation for the observed 
uniformity of the turbulence. 

5. Summary 

We have argued that energy to drive turbulence in 
the ISM of a galaxy can be extracted from the dif- 
ferential rotation either by way of a well-established 
MHD instability (Balbus & Hawley 1991), or directly 
from the Tr^ Maxwell stresses if the field is too strong 
to be formally unstable. In the bright inner disk, 
where stars are forming and the SNR is high, the 
classical picture of turbulence driven by stellar pro- 
cesses prevails, but where the stellar density is low, 
the energy deposition rate from MHD driven turbu- 
lence becomes important. 

A quite reasonable field strength (a few micro- 
gauss) is required for the turbulence to be Alfvenic, 
and its observed striking uniformity is, we argue, a 



simple consequence of flux freezing. But turbulence 
will decay without a constant energy source. We ar- 
gue that MHD instabilities are the preferred energy 
source in the outer H I layer because they are in- 
evitable and will maintain the characteristic turbu- 
lent velocity near the Alfven speed. Other sources of 
energy, SNe, infall, etc. are also available, and may 
be important locally, but do not provide a natural 
explanation for the magnitude and uniformity of the 
turbulence. 

By providing a mechanism for maintaining turbu- 
lence in the absence of mechanical energy input from 
stellar processes, we can claim some theoretical under- 
standing of the semi-empirical star formation thresh- 
old advanced by Kennicutt (1989). He argues that 
(significant) star formation in a galactic disk is trun- 
cated abruptly at the outer edge by the rising value 
of Q, assuming the observed constant velocity dis- 
persion of 6 km s^^. He shows that this simple crite- 
rion accounts for the extent of most star formation in 
many galaxies. Here, we are able to provide a possi- 
ble explanation for the constant velocity dispersion by 
linking it directly with an Alfven speed (cf. eq. [20]). 
We are unaware of any previous explanation for the 
maintenance of the observed turbulence required for 
Kennicutt's empirical rule. 

It must be stressed, however, that we have not ex- 
plained why the value a — 6 km s^^ should be so 
universal for disk galaxies. Linking it to the Alfven 
velocity implies a magnetic field of about 3 /iG - close 
to the field value inferred for many galaxies. Our 
achievement is to show that MHD instabilities can 
then maintain turbulent velocity fluctuations of the 
order of the the Alfven speed, but why this velocity 
should always be ~ 6 km s"^ is clearly related to why 
the magnetic field is always a few microgauss. This, 
we have yet to understand. 
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Table 1. Adopted properties in NGC 1058 



Radius 


5 kpc (100") 


10 kpc (200") 


H I column density 


6 X 10^" cm-^ 


4 X 10^" cm-2 


Total gas surface density Eg 


1.5 X 10~^ g cm~^ 


10"^ g cm~^ 


Full vertical thickness h 


400 pc 


400 pc 


Mean gas density p 


1.2 X 10-24 g cm-3 


0.8 X 10-24 g cm-3 


Circular velocity V^irc 


150 km s-i 


150 km s"^ 


Velocity dispersion in the gas a 


6 km 


6 km s"^ 


Epicyclic frequency k 


1.4 X 10-1^ s-i 


0.7 X 10-15 s-i 


Local stability parameter Qg . 


2.7 


2 


Jeans length Xrmcrit 


670 pc 


1,800 pc 
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